Symmetry breaking via orbital-dependent reconstruction of electronic structure in 

uniaxially strained NaFeAs 
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The superconductivity discovered in iron-pnictides is intimately related to a nematic ground 
state, where the C4 rotational symmetry is broken via the structural and magnetic transitions. We 
here study the nematicity in NaFeAs with the polarization dependent angle-resolved photoemission 
spectroscopy. A uniaxial strain was applied on the sample to overcome the twinning effect in the 
low temperature C2-symmetric state, and obtain a much simpler electronic structure than that of a 
twinned sample. We found the electronic structure undergoes an orbital-dependent reconstruction 
in the nematic state, primarily involving the d^^y- and dy^-dominated bands. These bands strongly 
hybridize with each other, inducing a band splitting, while the d^z-dominated bands only exhibit 
an energy shift without any reconstruction. These findings suggest that the development of orbital- 
dependent spin polarization is likely the dominant force to drive the nematicity, while the ferro- 
orbital ordering between d^jz and dyz orbitals can only play a minor role here. 



I. INTRODUCTION 

The parent compounds of iron-pnictides share a com- 
mon spin density wave or colHnear antiferromagnetic 
(CAF) state which is characterized by a ferromagnetic 
(FM) spin ahgnment along one direction in the two di- 
mensional rectangular lattice formed by iron sites and 
an antiferromagnetic (AFM) spin ahgnment along the 
other direction [InS] . Superconductivity in iron-pnictides 
emerges when the magnetic order is suppressed through 
chemical substitution or physical pressure [54^ • The 
mechanism causing the CAF order thus may be inti- 
mately related to the superconductivity. 

The development of the CAF order is always accom- 
panied by a lattice structure distortion, more specifically, 
a tetragonal-to-orthorhombic structural distortion. The 
transition temperature (Ts) at which the lattice distor- 
tion takes place either precedes or coincides with the CAF 
Neel transition temperature (T/v) The coincidence 

between the two transitions can be understood from a 
simple symmetry analysis. Besides breaking magnetic 
symmetry, the CAF state also breaks C4 lattice sym- 
metry of the tetragonal lattice. This lattice symmetry 
breaking due to the magnetic configuration can be sep- 
arately attributed to the development of a nematic elec- 
tronic order that only accounts for the rotation symme- 
try breaking but not the symmetry breaking of the mag- 
netic space. Based on the textbook Ginzburg-Laudau 
approach, the orthorhombic lattice distortion break the 
same C4 lattice symmetry. Therefore, the nematic elec- 
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tronic order and the orthorhombic lattice distortion are 
allowed by symmetry to be coupled with each other. 
Since the observed lattice distortion in all families of iron- 
pnictides is rather small, the electronic nematism could 
be dominant. 

So far, there are many experimental evidences support- 
ing the dominance of the electronic nematic behavior in 
iron-pnictides. Transport and neutron scattering mea- 
surements found a strong in-plane anisotropy of the re- 
sistivity and magnetic excitations [7l-[Toj. Angle-resolved 
photoemission spectroscopy (ARPES), scanning tunnel- 
ing microscopy, and quantum oscillation studies further 
revealed a complex electronic structure with C2 rota- 
tional symmetry pT| - [T5j . Moreover, the Fermi surface 
was suggested to be orbital polarized with almost pure 
orbital character [l^ in the CAF state. The en- 
ergy splitting of the dxz- and d^z-dominated bands was 
further observed far above T/v 12|. The observation of 
such a ferro-orbital ordering has led to a strong debate 
on the origin of the nematism: whether the electronic 
nematicity stems from the magnetic or orbital degree of 
freedom. Theoretically, some suggest that the nematicity 
is originated from spin fluctuations [l7l - [l9| , while others 
suggest the CAF state comes from a ferro-orbital order- 
ing |20l - [23| . Although following symmetry argument, the 
electronic nematicities from both origins must be coupled 
to each other, direct experimental evidence is needed to 
answer how strongly such an coupling is, and whether 
one is dominant over the other. 

Previous ARPES studies mainly focus on the so-called 
"122" series of iron pnictides, such as BaFe2--xCoj;As2 
and CaFe2As2. Although the mechanically detwinning 
methods were used to avoid the complexity due to the 
twinning effect in the CAF state, the observed electronic 
structures are still very complicated [ll| - [l3j . It is thus 
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difficult to obtain a compreliensive understanding of tlie 
complex electronic structure and its orbital character in 
the nematic state. In this paper, we report a systematic 
ARPES study on the uniaxially strained NaFeAs. The 
observed electronic structure is much simpler than those 
of the 122 series. Because of the polarization-sensitivity 
of orbitals in photoemission experiment, we could selec- 
tively probe the bands with different orbital characters. 
We found that the C4 rotational symmetry of the elec- 
tronic structure is broken via band reconstruction and 
the band reconstruction itself is also strongly orbital- 
dependent: the dxy- and d^z-dominated bands open a 
large hybridization gap at all the band crossings, while 
the dicz-dominated bands only exhibit an energy shift 
without any reconstruction. Moreover, the observed to- 
tal occupation of d^z is almost invariant. These results 
suggests that the ferro-orbital ordering between d^z and 
dyz can only play a very limited role in driving the elec- 
tronic nematicity. On the other hand, the orbital weight 
redistribution involving the dyz and d^y orbitals can only 
be explained through a strong coupling between spin and 
orbital degrees of freedom far above the Neel transition 
temperature, which suggests that the nematic spin fluc- 
tuations at high temperatures and their coupling with or- 
bital degree of freedom could be the driving force behind 
both the structural and Neel transitions in iron-pnictides. 



II. EXPERIMENTAL 

NaFeAs single crystals were synthesized by a self-flux 
method [23 | . Similar to other parent compounds of iron- 
pnictides, the high temperature phase (T > Ts) is in 
a tetragonal paramagnetic (PM) state, where each unit 
cell has two irons due to the inequivalent positions of 
arsenic ions in the FeAs layer [Fig. [T{a)]. Through the 
structural transition (oq > bo), there is an orthorhombic 
PM state (Ts > T > Ta?) with C2 rotational symmetry. 
The CAF magnetic order shows up when entering the 
orthorhombic CAF state (T < Tm) [IHi- The x- or k^- 
axes and y- or fc^-axes are defined along the AFM and 
FM directions, respectively. Previous ARPES study on 
NaFeAs observed a complex electronic structure in the 
orthorhombic CAF state due to twinned domains of the 
sample Taking the advantage of the mechanically 

detwinning method [8l.[25|. single domain with sufficiently 
large size was achieved in our experiment by applying 
uniaxial strain along Fe-Fe direction [Figs.[TIa) and[ljb)]. 
As a result, we could investigate the electronic structure 
of the orthorhombic CAF state without the complexity 
of twinning. 

The temperature-dependent resistivity of unstrained 
sample is consistent with the neutron scattering and 
transport measurements [1, H^. There are two inflex- 
ion points at 54 and 43 K, corresponding to the struc- 
tural transition temperature (T5) and the Neel transi- 
tion temperature (Tat), respectively [Fig. [IJc)]. Super- 
conductivity shows up at 20 K, which might be due to 
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FIG. 1: (Color online) (a) Cartoon of the lattice and 
spin structure in tetragonal PM, orthorhombic PM, and or- 
thorhombic CAF state. The x and y axes are defined along 
the iron-iron directions. The black arrows show the direction 
of the uniaxial pressure applied in the mechanical detwinning 
process, (b) The photograph of the device used to detwin the 
samples in our experiments, (c) and (d) The temperature 
dependent resistivity of unstrained and uniaxially strained 
NaFeAs, respectively. 



the doping effect of the deficient or excess sodium ions. 
When a uniaxial strain is applied, the in-plane anisotropy 
of the resistivity {pf,^ > pa^) emerges at about 75 K, 
much higher than Ts and Tn in the unstrained sample 
[Fig. [TJd)]. Similar behavior has been reported in uni- 
axially strained 122 series [8l-[To|. Recent neutron and 
X-ray scattering experiments have pointed out that the 
in-plane anisotropy of the resistivity in BaFe2As2 could 
be related to the lattice distortion [23, . The applied 
uniaxial strain could smear the structural transition and 
increase the onset temperature of the lattice distortion. 
Therefore, we attribute the 75 K observed here to be the 
structural transition temperature in uniaxially strained 
samples (Tg). 

ARPES studies were performed at the SIS beam- 
line of the Swiss Light Source (SLS). All the data 
were taken with a Scienta R4000 electron analyzer, 
the overall energy resolution was 15 meV and the an- 
gular resolution was 0.3 degree. The samples were 
cleaved in situ, and measured under ultra-high- vacuum of 
5xl0^^^torr. The polarization-sensitivity of orbitals in 
ARPES is a powerful tool to identify the orbital char- 
acters of band structure [1^, [sOj- The photoemission 
intensity is proportional to the matrix element of the 
photoemission process which can be described by 

J c>c [((/'jje • r|^|')|'^, where e is the unit vector of the 
electric field of the light, and cj)^ is the initial-state wave 
function. For high kinetic-energy photoelectrons, the 
final-state wave function cj)^ can be approximated by a 
plane-wave state e*'^ '" with k in the mirror plane. By 
rotating the polarization of incoming photons or the az- 
imuth angle of the sample, four experimental geometries 
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FIG. 2: (Color online) The experimental setup and matrix 
element calculations for three 3d orbitals. The mirror planes 
are defined by the analyzer slit and the sample surface normal. 
The two-dimensional plot of the Brillouin zone of NaFeAs is 
illustrated by red solid squares. Note that the photoemission 
cross-sections are amplified by a factor shown at the up-left 
corner of the panels. Thus, all the panels could be shown 
in the same color scales. There is a minor asymmetry in 
certain distributions caused by the out-of-plane component 
of the polarization. 



could be achieved named as p, s, p_r, and sn, respectively 
(Fig. [2])- The subscript "R" stands for the rotation of 
the sample. Since the low lying electronic structures at 
the Fermi energy (Ep) of iron-pnictides are commonly 
believed to be mostly composed of the d^z, dyz, and 
dxy orbitals (3l| . the matrix element distributions for 
these three orbitals were calculated in four geometries. 
As shown in Fig. [21 the matrix element distributions of 
the dxz and dyz orbitals exhibit a strong polarization de- 
pendence through out the first Brillouin zone, reflecting 
the opposite symmetry of these two orbitals. The ma- 
trix element distribution of d^y is suppressed along the 
direction parallel to the in-plane component of the polar- 
ization, due to its odd spatial symmetry with respect to 
the mirror planes. 



III. THE ELECTRONIC STRUCTURE IN THE 
TETRAGONAL PARAMAGNETIC STATE 

Based on previous studies [H,!!^, Fig. [3ja) sketches 
the Fermi surface of NaFeAs, which is consisted of two 
hole pockets (a and /3) around the zone center (F) and 
two electron pockets (ry and S) around the zone corner 
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FIG. 3: The sketch of the low-lying electronic structure of 
NaFeAs in the tetragonal PM state, (a) The sketch of the 
Fermi surface and its orbital character in the folded Brillouin 
zone for two iron ions per unit cell (solid lines) in the tetrag- 
onal PM state, (b) The observable Fermi surfaces in differ- 
ent experimental geometries, (c) The electronic structure of 
NaFeAs near Ep in the tetragonal PM state. 



(M) in the tetragonal PM state. The orbital charac- 
ters are dominated by the d^z, dyz and dxy orbitals [33| . 
The distribution of dxz is equivalent to the dyz orbital 
after rotating 90 degree, which is protected by the C4 
rotational symmetry in the tetragonal PM state. For 
simplicity, the continuous orbital evolution on the Fermi 
surface and the orbital mixing on the bands are not illus- 
trated here, which should be finite due to the hybridiza- 
tion of the bands with different orbitals. According to 
the calculation in section II, we could predict the observ- 
able Fermi surfaces in different geometries as shown in 
Fig.Efb). The polarization dependence of the Fermi sur- 
face is only sensitive to the in-plane component of the po- 
larization. The low-lying electronic structure of NaFeAs 
is illustrated in Fig. [3Jc). The bands dispersed along kx 
and ky directions are distinguished with the "x" and "j/" 
subscripts here and throughout the paper. As shown by 
previous study [3l|, the Px and ^x bands should con- 
nect with the band bottoms of rjy and 5y at the M point, 
respectively, which is guaranteed by the continuous evo- 
lution of bands in the unfolded Brillouin zone for one iron 
ion per unit cell. 

To further confirm the orbital character of the elec- 
tronic structure, the photoemission intensity maps ob- 
tained in p and s geometries are compared in Fig. Hlja) . 
Around F, the photoemission intensity distribution ro- 
tates 90 degree from the p to s geometry, which is con- 
sistent with the predicted distribution of the dxz and dyz 
orbitals on the a and /3 hole pockets. The electron pock- 
ets are observed either around My in the p geometry or 
around Mx in the s geometry, which is corresponding to 
the 5 electron pockets dominated by the dxy orbital. We 
did not observe the rj electron pocket here, which could be 
due to particular photon energy or experimental setup. 
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FIG. 4: (Color online) The polarization sensitivity of bands 
with different orbitals. (a) The photoemission intensity maps 
integrated over a 10 meV energy window at Ef taken at 80 K 
in the p and s geometries, (b) The photoemission intensity 
[I{k, cj)] taken at 80 K in p and s geometries along F-Ma; 
direction. The band dispersions are overlayed on the photoe- 
mission intensities with the solid lines, (c) The left and middle 
panels are the photoemission intensity maps integrated over 
a 10 meV energy window at 50 meV below Ep taken at 80 K 
in the p and s geometries, respectively. The observed cross- 
sections of bands are summarized in the right panel, (d) is 
same as panel b, but taken at 20 K. All the data are taken 
with 107 eV photons. 



since the outer rj electron pockets could be observed in 
previous studies [l^l or with 118 eV photons as shown 
later in Fig. [5j 

For the low-lying electronic structure, the bands dom- 
inated by dxz, dyz, and d^y could be selectively probed 
in different geometries [Fig. UJ^b)]. The ax band appears 
in the p geometry, while the fix band could be observed 
in the s geometry. The jx and Sx bands show weak pho- 
toemission intensities due to the small matrix element 
of the dxy orbital. These two bands could be observed 
more clearly in Figs. [5] and [6] We further investigated 
the photoemission intensity maps at 50 meV below Ep 
[Fig. life)]. Two ellipse-like features observed in p and s 
geometries are originated from the dxz and dyz orbital. 



respectively. They intertwined with each other forming 
the cross-sections of the a and /? hole pockets. The cross- 
sections of 7 is suppressed along the kx direction in the 
p geometry, which is consistent with the matrix element 
distribution of dxy When the sample is cooled down to 
20 K, the symmetry breaking of electronic structure is 
notable on the cross-sections of the a and /3 hole pockets 
[Fig. md)]. Both pockets extend along the AFM direc- 
tion. The dyz portions show more anisotropy than the 
dxz portions indicating an orbital-dependent reconstruc- 
tion of the electronic structure. 



IV. THE ELECTRONIC STRUCTURE OF THE 
NEMATIC STATE 

As shown in section III, although the simple matrix 
element calculation is not quantitative, there is a good 
qualitative agreement between the prediction and our ex- 
perimental results. Therefore, NaFeAs provides an ideal 
opportunity for us to study the temperature dependence 
of dxz , dyz , and dxy separately along both AFM and FM 
directions. We will focus on the temperature dependence 
of the bands that make the electron pockets first. 



The temperature dependence of the rj and 5 
bands 



There are two electron pockets r] and S around the zone 
corner. The rjx bands could be clearly identified crossing 
Ep along the kx direction in the p geometry [Fig. ^a.)]. 
Their Fermi crossings move outward continuously as the 
temperature decreases [Fig. [5jb)]. The momentum shift 
of the rjx Fermi crossing is about 0.02 7r/a from 60 to 
20 K. For the 6 electron pocket, we took the data along 
the momentum cut around My in the p geometry in or- 
der to enhance its photoemission intensity. As shown 
in Figs. [5jc)-[5l[e), when the temperature decreases, the 
MDCs peak intensity of Sx suppresses near Ep promptly, 
indicating an energy gap opening for the Sx band along 
the kx direction. Note that, the band dispersions along 
the kx direction should be equivalent for Mx and My, 
which is protected by the translation symmetry. 

Along the ky direction, the photoemission intensity of 
i]y is too weak to be distinguished [Fig. ^i)]- As will 
be discussed later, this is because the riy band shifts up 
above Ep. The tiny MDCs peaks observed in Fig. ^g) 
might be due to the minor twinned domains of the sam- 
ple. For the Sy portions of the inner electron pocket, 
weak temperature dependence was observed [Figs. [Sfh) 
and Elji)] . The 5y Fermi crossings slightly shrink with a 
momentum shift of about 0.02 7r/a from 60 to 20 K. 

The in-plane anisotropy of the electron pockets pre- 
serves at the temperature as high as 60 K, which is 
far above the structural and magnetic transition tem- 
peratures in unstrained sample. This suggests that the 
anisotropic electronic structure observed here might be 
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FIG. 5: (Color online) The temperature dependent photoemission data of the rj and S bands around the zone corner, (a) The 
I{k,uj) taken at 60 and 20 K in the p geometry along Cut #1. (b) The temperature dependence of the momentum distribution 
curves (MDCs) at Ep along Cut #1. The momentum position of Cut #1 is shown by a black arrow in the inset of panel b. 
(c) The I{k,Lo) and its second derivative image with respect to energy [9^/(fc, a;)/i9a;'^] taken at 60 and 20 K in the p geometry 
along Cut #2. (d) The temperature dependence of the MDCs at Ef along Cut #2. The momentum position of Cut #2 is 
shown by a black arrow in the inset of panel d. (e) The MDCs selected between -36 meV to 3 meV for the data in panel a. 
Each MDCs has been individually normalized by its integrated weight to highlight weak features, (f) and (g) are the same as 
panel a and b, but taken in the pr geometry along Cut #3. (h) and (i) are the same as panel c and d, but taken in the pr 
geometry along Cut #4. (j) The sketch of the electron pockets in the tetragonal PM and orthorhombic PM states. All the 
data are taken with 118 eV photons. The color scale of I{k,ui) and d'^ I{k,ui) / dui^ are shown in the bottom-right inset. These 
color scales and notations are used throughout the paper. The band dispersions are overlayed on the photoemission intensities 
with the solid lines, which are determined with the help of the second derivative images and MDCs. 



closely related to the strain and thus lattice distortion, 
whose onset temperature {T'g) increases due to the uni- 
axial strain (More evidences in Fig. |5]). Therefore, in the 
orthorhombic PM state above Tjv, the electron pockets 
reconstruct significantly as shown in Fig. [SJj). Two por- 
tions vanish in the the orthorhombic PM state and the 
remaining electron pocket extends along the kx direction. 

B. The temperature dependence of the a, P, and 7 
bands 

To investigate the anisotropy of the /? and 7 bands, 
temperature dependent data taken in the s geometry 
along the kx direction are shown in Fig. [SJa). Remark- 
able band reconstruction was observed. Near F, the /3x 
and "fx bands open a large hybridization gap at 20 K. 



As a result, the Fermi crossings of /3x are pushed out- 
ward from F, and the band top of 7 shifts downward 
to about 40 meV below Ep [Figs. EJa), Eljb), andi^d)]. 
The f3x band disperses towards Ep around M forming a 
hole-like band. It hybridizes with the Sx electron pocket 
and opens an energy gap at Ep upon lowing the tem- 
perature. This is consistent with the gap opening of Sx 
observed in Figs.[5][c)-[5{e). Moreover, the band shift and 
hybridization push the Px band above Ep around M. 
Considering the fact that the /3x band should connect 
with the band bottom of rjy [Fig. [Sfc)], we could further 
confirm the observation in Fig. [5l^f) that the r]y portions 
of the outer electron pocket shift up above Ep in the 
orthorhombic PM state. Note that, the long-range mag- 
netic order could halve the Brillouin zone by folding the 
bands between F and M. We thus attributed ry^ near F 
together with /3'^ and a'x around M to be the folded bands 
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FIG. 6: (Color online) The temperature dependent photoemission data of the /3 and 7 bands, (a) The I{k, lj) and d^I{k, uj)/duj^ 
taken at 60 and 20 K in the s geometry along Cut #1. (b) and (c) The Energy Distribution Curves (EDCs) selected from 
the data in panel a around F and M, respectively, (d) The temperature dependence of the MDCs at Ef along Cut #1. The 
momentum position of Cut #1 is shown by a black arrow in the inset of panel d. (e) and (f) are the same as panel a and d, 
but taken in the sr geometry along Cut #2. (g) The sketch of the /3 hole pocket in the tetragonal PM and orthorhombic PM 
states. All the data are taken with 118 eV photons. 



in the magnetic state. 

In contrast to the complex behavior along k^, the 
/3 band does not reconstruct along the ky direction 
[Fig. mje)]. The Fermi crossings of /3y are almost tem- 
perature independent as shown in Fig. [S^f). Around the 
M point, the (iy band shifts downward to about -50 meV 
at the lowest temperature. Since (3y should connect with 
the band bottom of rj^ [Fig. Efd)], such an energy shift 
of f3y is consistent with the enlargement of rjx observed 
in Fig. [Slja). We summarized the reconstruction of /3 in 
Fig. IHg). The P hole pocket extend along the AFM di- 
rection in the orthorhombic PM state. 

The reconstruction of the inner a hole pocket is shown 
in Fig. [71 Along the ky direction, the ay band split into 
two parts indicating possible hybridization with other 
bands [Figs. [71[a) and[7l^c)]. The hybridization with S'y 
could be first excluded. According to the temperature 
dependent data, while band shift of the lower part of ay 
is notable at 60 K [Fig. [ZUd)], the folding of Sy occurs 
at a much lower temperature of about 49 K [Fig. El^b)]. 
Moreover, no hybridization was found for the 6y band 
in Fig. EJh). Therefore, the ay band should hybridize 



strongly with the 7^ band in the orthorhombic PM state, 
as shown in the inset of Fig. El^c). In that large 
proportion of the dyz orbital on ay would mixed into "fy 
due to such hybridization. Consequently, we note that 
the photoemission intensity of 7^ is suppressed around 
the r point at 20 K in Fig. ^e), since the dyz orbital 
could not be observed in the sr geometry. Along the 
kx direction, the bands show weak temperature de- 
pendence. Their Fermi crossings slightly move outward 
about 0.02 7r/a from 60 to 20 K [Figs.[7i;e)-17i;g)]. There- 
fore, similar to the /3 hole pocket, the a hole pocket ex- 
tends along the AFM direction in orthorhombic PM state 
as shown in Fig. [D^h). 

Our previous APRES study on unstrained NaFeAs 
has found that the short range magnetic order emerges 
just below the structural transition temperature of 54 K, 
which causes a band folding [23|. However, when the 
structural transition is enhanced to 75 K here by the 
uniaxial strain, we still observe the folding of Sy bands 
below ^49 K [Fig.[7I^b)]. Therefore, the short range mag- 
netic order observed above T/v in previous study should 
be related to the fluctuation of the long range magnetic 
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FIG. 7: (Color online) The temperature dependent photoemission data of the a band around the zone center, (a) The I(k,bj) 
and i9^/(fc, a;)/i9a;^ taken at 60 and 20 K in the p_R geometry along Cut #1. (b) The temperature dependence of the MDCs 
at Ef along Cut #1. The momentum position of Cut #1 is shown by a black arrow in the inset of panel b. (c) The MDCs 
selected between -72 meV to 6 meV for the data in panel a. Each MDCs has been individually normalized by its integrated 
weight to highlight weak features, (d) The temperature dependence of the MDCs at 50 meV below Ef along Cut #1. (e), (f), 
and (g) are the same as panel a, b, and c, but taken in the p geometry along Cut #2. (h) The sketch of the q hole pocket in 
the tetragonal PM and orthorhombic PM states. All the data are taken with 118 eV photons. 



order, instead of a direct correlation with the structural 
transition. 



V. DISCUSSION 

To further investigate the orbital-dependent behaviors 
of the reconstruction in the orthorhombic PM state, the 
bands dominated by the and dxy orbitals are high- 
lighted in Fig. [8ja) . The most remarkable feature found 
here is the hybridization gap opening at almost all band 
crossings. As a result, the bands split into the upper and 
lower bands around F and open an energy gap at Ep 
around M. Such reconstruction could redistribute the or- 
bital weight of dyz and dxy, inducing a partial gap at Ep 
[Fig. mjb)]. On the contrary, the da;z-dominated bands 
exhibit no remarkable reconstruction, except an energy 
shift around M [Fig.[8jc)]. The orbital weight redistribu- 
tion is negligible for the d^z orbital. Since the energy shift 
of the /3 band is pronounced near the M point, we use this 
energy shift as a parameter to describe the reconstruction 



of the electronic structure. If the C4 rotational symme- 
try of the electronic structure is preserved, the /3 band 
dominated by the dxz and dyz orbitals should degenerate 
along the kx and ky directions. As shown in Figs. IH^d)- 
[5{f), when the temperature decreases, the energy split- 
ting between fix and /3y bands enlarge linearly and almost 
saturate at T/y. The onset temperature of such splitting 
could be estimated to be around 75 K, where the in-plane 
anisotropy of resistivity emerges This result is con- 

sistent with previous studies on the "122" series of iron 
pnictides 12], and suggests that the significant orbital 
weight redistribution occurs below the structural transi- 
tion temperature. Moreover, the smooth evolution of the 
electronic structure reconstruction into the CAF state in- 
dicates the both the magnetic and structural transitions 
share the same driving force (23 |. 

As proposed by many theories, the orbital weight redis- 
tribution is important to understand the nematic tran- 
sitions in iron-pnictides [211-23, 34 37]. Our results thus 
provide several important quantitative tests for these the- 
oretical descriptions. We could first check the inequiva- 
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FIG. 8: (Color online) The orbital-dependent reconstruction 
of the electronic structure, (a) The low-lying electronic struc- 
ture in the tetragonal PM and orthorhombic PM state re- 
spectively, where only the dyz- and d^iy-dominated bands are 
highlighted. The band dispersions above Ef are shown con- 
sidering the continuous dispersion of the bands with certain 
orbitals. (b) Cartoon shows the orbital weight distributions 
near Ep in the tetragonal PM and orthorhombic PM states, 
(c) The low-lying electronic structure in the tetragonal PM 
and orthorhombic PM state respectively, where only the d^z- 
dominated bands are highlighted, (d) The temperature de- 
pendence of the EDCs at k\ as indicated by the gray line in 
panel a. (e) The temperature dependence of the EDCs at 
as indicated by the gray line in panel c. (f ) The peak positions 
of the /3a; and fiy bands as functions of temperature. 



and dyz orbitals, which 



lent occupation between the dx 
is considered to be important for the formation of ferro- 
orbital ordering [2l| - |23j . The occupation change of the 
dxz orbital could be calculated from the size increase of 
the r\x electron pocket subtracted by the size increase of 
the ax hole pocket. The momentum shift of the r\x Fermi 
crossings is about 0.02 7r/a from 60 K to 20 K. Consider- 
ing the band shift above 60 K [Fig.[51[f)], this momentum 
shift might be doubled to 0.04 7r/a corresponding to a 
maximum occupation change of about 0.04 electron per 
iron. Moreover, the Fermi momenta variations are almost 
the same for the enlargement of the r\x electron pocket 
and the ax hole pocket. Therefore, the occupation change 
of the dxz orbital should be infinitesimal, which suggests 
that the ferro-orbital ordering only involving dxz and dyz 

ts a minor role in driving the nematic transitions here 
Note that, although the occupation of dxz is almost 
invariant, the orbital polarization could be strong only at 
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FIG. 9: (Color online) The temperature dependence of the 
Fermi surface morphology, (a) The photoemission intensity 
maps in the orthorhombic CAP state taken around F and M 
in the p geometry, (b) is the same as panel a, but taken in 
the s geometry. The intensity was integrated over a 10 meV 
energy window at Ef- (c) The Fermi surface morphology in 
the orthorhombic CAF state, (d) The band structure along 
kx and ky directions in the orthorhombic CAF state, (e) The 
summary of Fermi surfaces in the tetragonal PM, orthorhom- 
bic PM, and orthorhombic CAF state respectively. 



the low energy. Since the orbital weight of dyz and dxy is 
partially gaped away from Ep, the Fermi surface could 
exhibit more dxz orbital character in orthodromic state 

MM- 

More intriguingly, instead of the pure ferro-orbital or- 
dering between dxz and dyz orbitals, the reconstruction 
of the electronic structure involves mostly the dyz and 
dxy orbitals. These two orbitals are found numerically to 
contribute much larger magnetic moment than dxz in the 



magnetic state 



Therefore, the orbital weight 



redistribution observed here below Tg could be related to 



the spin polarization of the dyz and dxy orbitals. 



In this 

scenario, the strong coupling between spin fluctuation 
and orbital degree of freedom could induce an orbital- 
dependent spin polarization, which further reinforces the 
spin fluctuation and the anisotropy of electronic structure 
far above the Neel transition |34l437| . Such a positive 
feedback process can be accelerated by the an uniaxial 
strain, which helps to explain the dramatic enhancement 
of Ts by the rather small strain applied in our experiment 
Q. Consistently, the enhancement of spin fluctuation 
could be observed under the structural transition js^. 
Our results therefore suggest that the nematic spin fluc- 
tuations at high temperatures and their coupling with or- 
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bital degree of freedom could be the driving force behind 
both the structural and Neel transitions in iron-pnictides. 

In the magnetic state, the electronic structure is folded 
between F and M due to the long-range magnetic order. 
According to the photoemission intensity maps in s and p 
geometries [Figs. |9l[a) and[9l[b)], the Fermi surface mor- 
phology and band structure in the magnetic state are 
shown in Figs. [Sljc) andlHlJd). The electron pocket folds 
from M to F and then hybridizes with the /3 hole pocket. 
As a result, four small pockets are formed. There are 
two bright spots along AFM direction dominated by the 
dyz and dxz orbitals. Since they are formed by two bands 
with opposite symmetries [Figs.|9l[c) and|9l[d)], these two 
small pockets are intensively discussed as Dirac points in 
122 series [Sjll^Eol. The other two pockets are formed 
by the Sy and /3y bands along the FM direction [Figs.jSJc) 

We summarized the reconstruction of Fermi surface 
in Fig. inje). The C4 rotational symmetry of the Fermi 
surface is broken under the structural transition. The 
anisotropic Fermi surface observed in the the orthorhom- 
bic PM state could be responsible for the in-p lane 
anisotropic transport behaviors above Tjv [sl-fiol |41|. All 
the pockets show anisotropic shape with an extension 
along kx direction. Especially, the transport proper- 
ties contributed by the electron pockets could be highly 
anisotropic, since two portions of electron pockets are 
gaped away from Ep in the orthorhombic PM state. 
When entering the orthorhombic CAF state, the Fermi 
pockets are folded between F and M. The Fermi surface 
is thus characterized by four distinct pockets originated 
from the hybridization between the hole and electron 
pockets. These features qual itativel y ag ree with previous 
studies in the 122-series [ill , [l^ [H, H^, which highlight 
the unified nature of the nematic magnetic state in iron- 
pnictides. Furthermore, the Fermi crossings of /?, Sy, and 
r]x are sensitive to the structures, kz, and doping lev- 
els. Therefore, the topology of these four pockets could 
be variable in different compounds. For example, the 
pockets distributed along FM direction are hole-like in 
NaFeAs, but electron-like in 122 series [li, I2„ 15., ^] . In 
that case, the topology change of these distinct features 
could be responsible for the non-monotonic doping de- 
pendence of the in-plane resistivity anisotropy observed 



in 122 series [i-[ll|42. 



VI. SUMMARY 

To summarize, our polarization-dependent ARPES 
study on uniaxially strained NaFeAs has clearly demon- 
strate the details of the electronic structure reconstruc- 
tion in the nematic state. We found that the rotational 
symmetry breaking of electronic structure is closely re- 
lated to the in-plane anisotropy of the resistivity, whose 
onset temperature is far above the Neel transition tem- 
perature. More intriguingly, such a symmetry break- 
ing originates form the reconstruction of the electronic 
structure, and exhibits orbital-dependent behaviors. The 
dyz- and d^jy-dominated bands split near F and open an 
energy gap around Af, while the daj^-dominated bands 
only exhibit an energy shift without any reconstruction. 
Our observations exclude the theories based on the ferro- 
orbital ordering only involving the d^z and dyz orbitals, 
since the occupation of d^z is almost temperature inde- 
pendent. Instead, the strong coupling between the ne- 
matic spin fluctuations and the orbital degree of freedom 
are suggested to drive the nematicity. Such a coupling 
causes an orbital-dependent spin polarization, and thus 
the development of finite magnetic moments of the dyz 
and dxy orbitals as manifested in their weight redistribu- 
tions. Our findings provide the experimental foundation 
for a microscopic understanding of the nematic state of 
iron-pnictides. 
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